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ABSTRACT 


The magnetotelluric (MT) method of sounding is 
reviewed and the interpretative techniques including 
both numerical modelling and tensor impedance concepts 
are presented in a formulation suitable for the study 
of the problem of noise in magnetotelluric signals. 

A statistical treatment of the noise problem 
has resulted in the development of a cyclical process 
rater produces a pronounced enhancement of the signal 
to noise ratio over the complete spectral range in the 
Majority of processed data samples. In those cases 
where the data was not so improved a definite confidence 
limit could be placed on that portion of the spectral 
range which was acceptable. 

This technique was applied to magnetotelluric 
soundings made at 14 sites around the Black Hills of 
South Dakota in the frequency band 10a ComLOrMHZe “rhe 
soundings cover an area of approximately 30,000 square 
Km and correspond to a depth range of 0.1 to 50 km. 

A three dimensional earth model of the estimated 
resistivity distribution is derived from the sounding 
results, using the generalized impedance concepts. An 
anomaly of high conductivity at depths of 3 to 10 Km is 
found in the central region of the Hills, under the 


cover of highly resistive granitic material. 
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CHAPTER i 


INTRODUCTION 


1-A. Historical Review 


The process of determining the subsurface resis- 
tivity structure of the earth, by measurements made at 
the surface, using the naturally occurring electro- 
magnetic fields as a source, was originally proposed 
by Tikhonov (1950). Measurements were carried out 
both by Tikhonov (1950), Kato and Kikuchi (1950), and 
Rakucakeve(1950 C951) ~“Inta) definitive work?) Cagniard 
(1953) presented a graphical method for the interpre- 
tation of results in the case of a horizontally layered 
earth and introduced the expression "Magnetotelluric 
Method" to designate the technique. 

As originally envisaged by Cagniard, the source 
was infinitely extended and thus a plane wave impinged 
on the horizontally stratified earth. Wait (1954), 
and later Price (1964), investigated the effect of 
finite source dimension on the measured impedances, 
and argued that if the source is finite in extent or 
possesses horizontal gradients, complex angles of 
incidence must be considered. However, Madden and 
Nelson (1964) showed that for reasonable earth con- 


ductivity models, the plane wave assumption appears to 
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be valid in the frequency range of interest for the 
MT method. 

A more important consideration than that of 
source dimensions, may be that a one-dimensional 
model is inadequate for a real earth where lateral 
inhomogeneity and anisotropy in the conductivity 
are present, and the electric and magnetic fields 
may, in general, no longer be orthogonal. d'Erceville 
and Kunetz (1962) and Rankin (1962) gave analytic 
solutions for a vertical fault and vertical dyke 
respectively. Chetaev (1960), Cantwell (1960), 
Bostick and Smith (1962), Mann (1965), O'Brien and 
Morrison (1967), Praus and Petr (1969), Rankin and 
Reddy (1969, 1970), and others have studied various 
types of anisotropy in conductivity by means of 
impedance tensor technique, and have proposed methods 
COmnTiICmEnewpLInctpal directions, OL Conductivity 
anisotropy. Madden and Nelson (1964), Swift (1967), 
Morrison et al. (1968), Sims and Bostick (1969), Word 
etpal. = (1970) Vozort (19/2), -andwothers have presented 
technique to compute tensor impedance elements; this 
tensor matrix relates the electro-magnetic fields 
for two and three dimensional structures. Analog 
model studies have been carried out by Rankin et al. 


(1965), Dosso (1966), Takaes (1969), and others. 
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A considerable number of experimental results 
have been reported from various locations in North 
America, Europe, and Russia (Srivastava et al. 1963, 
Whitham and Anderson 1966, Berdichevsky 1966, Vozoff 
1969, 1972, Peeples and Rankin 1973, and others). 
While the results have been significant, many of the 
questions and problems involved, particularly with 


respect to interpretive techniques, still remain. 


I-B, Outline of the Thesis 


Chapter two will present the theory of the MT 
methods. 

Chapter three will describe the recording system 
and the techniques of data analysis. 

Chapter four will discuss the noise effect, and 
present a technique for improving the signal to noise 
seSli ea ele 

Chapter five will present the sounding results, 
geological background, and model interpretation for 


the Black Hills in South Dakota. 
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THEORETICAL DISCUSSION 


II-A. Maxwell's Equation for Plane Waves in a Uniform 


(homogeneous and isotropic) Medium 


Consider the model in the cartesian coordinate 
system represented by xyz-axes as shown in Fig. (2-1). 
The z-axis is verticallydownward and the earth of 
conductivity 9(x,y,z) occupies the half space z > 0 
with the surface at z = 0. The source is in free 
space z<0Q. Maxwell's Equations in the MKS system 
Jwt 


for e time dependence are 


vx B= - B= -jupH (2-1a) 
at bie ah : 
> > > 
Vx H= J = (o + jwe)E (2-1b) 
> 
V- E= q/e (2-1c) 
> 
V-« H= 0 
where w = 27f is the angular frequency; gq is free charge 


density; and o,u,€ are the conductivity, permeability, 
and dielectric constant respectively, of the media. 
Except for anomalies which will not be discussed in 
this thesis, u = ey, and € = Eo: For a uniform half 


space, the wave equation can be obtained from (2-1a,b) 


by taking the Curl. 
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Vx Vx A = -jwulo + jwe)A (2-2) 


> > > 
where A represents a field vector E or H. Eqn. (2-2) 


describes the field for all regions. 

Except at discontinuities or in the source 
region, the free charge density is zero and hence 
within the medium V-E = 0 and Eq. (2-2) becomes the 


Helmholtz Equation 


v-R = jJwulo + jwe)A , (2-3) 


In a homogeneous region where o is constant, Eqn. (2-3) 


is satisfied by a plane wave solution of the form 


> > 
A=A af cee 
fo) 
where 
r= sie 4 fy + kz is the position vector, and 
k = ik,+ jk, + kk, is the propagation vector, and 
ae 
Kee Ki = Vjup (ort ajwe) Ce) 


The frequency of interest in this work will 


range from 10°? — 10 Hz, and the conductivity of 
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the earth from 107 
the earth, and henceforth the propagation constant 


in the earth becomes approximately 
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In the air where o = 0, the propagation constant is 
K, = jwvue d (2-5b) 
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When the field propagates in the earth, the real 
component of Ky produces an attenuation effect on the 
field. The attenuation factor for a given wave in a 


medium with conductivity o is 


Re[K, ] = /wud/2 or TU 


The propagation distance or penetration depth 6, within 
which the strength of the field is attenuated by a 


Pacloneos soe, is defined as the skin depth, where 
6 = 1/Re[K, ] = myeer (2-6) 


The wave length of the field in this conductive medium 


is related to the skin depth by 
A = 270 . 


The field measured at the surface due to reflect- 
ing surfaces at a depth equivalent to approximately one 


half-wave length for example, will be attenuated by a 
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Eactor a cae In addition the secondary fields are 


reduced by the finite reflective coefficient of the 
anomaly and possibly out of phase relations. 

The velocity of phase propagation of the wave 
WiLtNanwal Meal UuMmeOfbeCONGUCTIVILY 40, 1S GAs elnus. for 
a conductivity o = 0 ke maOoseand cma ul0NeeH2; athe 
velocity is 

v= + fag hes 107 m/sec . 

As a consequence a plane wave incident at an oblique 


angle will propagate perpendicular to the surface as 


though it were at normal incidence. 


II-C. Plane Wave Incident on a Uniform Earth 


In the case of a uniform earth where o is cons- 


tant, the E and H-fields are related in Eqn. (2-la) 


Considering the E-field polarized in the x-direction 
and the H-field in the y-direction, the above relation 


can be rewritten as 


OEY 
Bre bel eee or okies 
and 
E. = Yjwup H (2a) 


Hy s 
( : 7 1 
hs ; eam Re, 
on? a One 
i Vv et 
a ae ai ae = oe 
W3EB 404 Y vrabacse@ Se marae a ec 
ao 20 iokiisag avrdusd yor, S32 ne oi ei Beanibex 
tteloy eaeay to; I8O Ys otha | aR eistions 
Lutgepaiord Seat to we ivetan ae = 
. ‘ 
gi ,O YRLVEISeb rote, Te pau ty anise vi - 
, ot i " 
. Bris act Uf =.6 Ne ivivoulacd & 
as, al 
1 as 
$ 7 a 
a a Lm i 
- > a \ 
aris Bi. 2 HOD SD til BW. said rege soe Late TE LS Mee) 
7 = ; {hy 
5 Lome ie PSOE. b. ese 
SO>aGi) 1 t Gn oe to alee 2B ow: bolle 
ij od HEY OG he 3 for VE Ae 
> — - ~ a a> = Fe ieee eR np eng A a ee ey 
iss!) muotivw, 5 26 Shes eae wee 
pe po laries bre 2 aad 
re oe Aerie 
=) . Bes = Hs a 
nots pace x ar? ae ae a wee ook, beara pire 


/~ ’ 
oe 
- 7 


poner chr itte eid zd oeaAbE oa s 


hs 
S6 


re i 
a o in 


es 


where 90 = 1/o is the resistivity of the medium. 


The impedance Z is defined 


N 
Wt 


E ‘ 
w= Yjoup = Youp 2" (2-8) 


Thus the phase angle between the E and H-fields for 
a uniform half space is ¢ = 71/4. The resistivity in 
this case of the uniform earth obtained from Egn. (2-7) 


Oreale- Seals 


° (2-9) 


II-D. Apparent Resistivity of Layered Earth 


For a uniform earth as discussed in the previous 
section, the true resistivity of the medium can be 
directly deduced from the measured fields. For a 
more general case, however, the true resistivity of 
the medium is a more complex function of the measured 
fields. 

Ce) single ‘layered earth: -This model consists of a 
uniform layer of thickness h and resistivity Py, over- 
lying a half space with resistivity Por as shown in 
Fig. (2-2a). Choosing the source E-field polarized in 
the x-direction, the plane wave equation can be written 


(from Eqn. (2=3)))— as 
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A solution has the form 
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in each layer. The coefficient A; is associated with 
a poSitive going (downward) transmitted wave, and B. 
with the upward reflected wave. Zs4 is the depth 

to the upper surface of the i-th layer. Thus in 


medium l, where Z. =z = 0 
hdl 


HAZ) = A 


Ther Sabge- 9 
E, (2) = K,p, [A,e Bie ] 


while in medium 2 where Zee ee h (the thickness of 


the first layer) 


-K, (2-h) 
e 


HAZ) A, 


~K, (z-h) 
E, (z) = K,p5Ae 
where Bo = 0 since there is no upward reflected wave. 
From the boundary condition that the tangential com- 


ponents of the field are equal at z =h, 
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Thus, the surface impedance is obtained 
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where a= e ° 
The apparent resistivity is defined by 
A, -B, }2 
ay ah Zi ll i s 
P. = mn Raee Timi ae - (2=1'1)) 


This is the apparent resistivity of a single-layer 
earth, which differs from the true resistivity Py by 
amLaccor (A, -BDAA,+B)1°. This) factoreis a f£unction 
of the earth parameters, i.e., the resistivities of 
the two differently conductive regions and the thick- 
ness of the overlying layer. 

(2) n-layered earth: A typical n-layered earth 
model is as shown in Fig. (2-2b). The solutions for 


the m-th layer have the form 
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Applying the boundary conditions at z = Za? the trans- 


mission and reflection coefficients can be expressed 


by 


iL Rm 
A = xLU+tyA Ly “ (1-y)B ile 


ih am 
B= gl (-YA4, + (-y)B,,,Je 


where y = (K See) ewe lee Zn-] 18 the 


sen ens m 


thickness of the m-th layer. The impedance of this 


m-th layer is defined 


E Kem tani (hs >of ) ato 
Ve ps m m mm 
m H mm kK 


eee sa 5 (2212) 
y,m mm m+1 


tanh ee 


The surface impedance is 


sie B (0) sip A, - seas . Ki Py tanh (K,P,) + Zo 
fe) H,, (0) eae A, +B, ay KP, + 2, tanh K,P,) 
(2-13a) 


and the apparent resistivity is 


A,-B 
2 Tima 
Lio ey (2-13b) 

1|A, +3, 


' 2 
where Aj Bi and ew iBereuEn [A,-B,/A,+B,] are functions 


of the resistivities and thickness of all the layers. 
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II-E. Tensor Impedance Analysis 


In the previous section, the case of a uniform or 
piecewise uniform medium was discussed in which the 
field component variations were a function of one 
variable only and the resistivity was a scalar. In 
a more general case both the components of the fields 
and the resistivity are functions of 3 spatial variables. 
In this case one defines an impedance tensor [Z] such 


that 


> > 
[BE] = [2] [H] (2-14) 


where [E] and [H] are the electric and magnetic field 
vectors written as column matrices and [Z] is in general 
a 3 x3 impedance matrix relating each of the components 
of E to 3 components of H. The components of [Z] are 
point functions of earth parameters together with the 
appropriate frequency dependence. While analytic solu- 
tions are possible only in very special cases, it is 
of interest to develop general expressions for the 
electromagnetic fields in terms of the earth parameters. 
Assuming the earth behaves like a linear system, 
the secondary field and in turn the total field is 
linearly related to the primary magnetic field, at de 


At the surface, x,y,z = 0 


H 2 ats 
H(0)] = (HL (0), HA (0), 0) (2-15) 
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Thus, the total H-field within the earth can then 


be related to the primary field at the surface 
> > 
ei PN) SP roa a Eas ibeh hal (2-16) 


where [o(r)] is generally a 3x3 matrix. However, since 
Hz (9) = 0, only the 3x2 part is present. The elements 
of [o¢(r)] are functions of the transmission and reflec- 
ElLonecoerretcients, which are in turn, of the position r 
and other earth parameters. 

The E-field at position r can be obtained from 
a generalized form of Ohms law, appropriate to an aniso- 


tropic earth, and Maxwell's equations, 


[E(r)] = [o(r)] (3 (r)] = [p (x) IV x H(x)] (2-17) 


where op(r) is a 3X3 matrix, its elements representing 
the tensor coefficients of the resistivity of the 
earth; and J is the current density. At the surface, 
the components of J in terms of rae can be obtained 
by substituting H(r) from Eqn. (2-16) into J.(r) = V x Hr) 


and evaluating at the surface 
0 Ta ky is ine ade ote games is Demi aay yet eae 


ae 2- 8b 
tly (2—-18b) 
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where _$.. = (oe erat )¥) : 


ch ta inal r=0 

Since no conduction current flows 
region, the continuity condition at the 
interface requires that in the earth at 
surface, J, vanishes. Thus the E-field 


can be expressed in terms of ete by 


E, (0) P4495, (0) + Pyoty'o) 


Hoy 69) (2-53.G) 


inesche air 
air-earth 
or near the 


at the surface 


Sey Paes) Sel P15 7945 TH. (0) 


+ 10144939 - 220) * Ppa loor27 xb32) Hpy (0) (2-19a) 


Sus Po 7Fy (0) aS Po 24, (9) 


= [p97 49317 2921) 7 ©22 ‘24117 x31 


+ 


E,(0) = 9315, (0) + 9359, (0) 


Lpoy (yo327 2922) * P22 '2o12 ~xb32) I Bpy (9) 


1H, (0) 


(2—19b) 


= [031 (y$31 7 2921) * 932 (221 ~ 4932) Hp, ©) 


Ploary?somize22) | Pao 291d x32) apy) 


(2-19c) 
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Eqns. (2-19) define the primary impedance relating the 


Primary H-field to the E-field at the surface, by 


PC) a [2 1TH, (0)] (2-20) 


where [2] is a 2x 3 matrix as indicated in Eqns. (2-19). 

The surface impedance, which defines the relation- 
ship between the total H-field and the E-field at the 
surface, can be determined as follows. From Eqn. (2-16), 
the inverse matrix of [4(0)] can be generated and 


evaluated at the surface, 


i 


[H, (0) ] =e (0) | ee alH (0) |] (2-215) 


where [6(0)]7+ is defined by [4(0)]~*[(0)] = [I]. 
SubsSLIcutLingeton (2-21) tom 2—-20) ,athe total H-field 


and the E-field at the surface are related by 


[E(0)] = [2,1 iakiips HG) (2-22) 


Thus, the surface impedance function for the three 


dimensional earth is 


Uh 


ra) = [2] [o(0)] (2-23) 


where [Z] is a 3x 3 matrix. An explicit form of E-H 
fields relation can then be expressed in terms of the 


elements of [Z]. 
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BEY TAR ali) 9 Dae! S/ AE Pal 9 get ey AE | (2-24a) 


H (2-24b) 


E = 24 H.. ae gs Ht “Ape ea | (2-24c) 


IIi-F, Examples of the Impedance Functions 


(a) Unpeornmwarths ae o(X,y,2) = oh 1/p.,)- 


In a uniform half-space the total H-field is the 
transmitted primary H-field. If the surface value of 


the field is denoted by 


> 
H(0) = (Ho, Ho, 0) , 


then the total field at z>0 is 


Be ae = —KZ —KZ 
H(z) = Ne) = (Hoe ; Bho a EO 
for which, 
Cue 0 
[o(r)] = |0 ae (a5) 
0 0 


where Kev WHO, ° 


Thus, at the surface, the inverse of [9(0)] is 


1 0 0 
iy Gh alae | | : 
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The resistivity of a uniform earth can be denoted by 


0 
[pop] = |0 p 
0 0 e) 
and so that 
0 oK 
Z = |-oK 0 ‘ 
Med le 
0 0 


Therefore, the surface impedance tensor is, 


0 oK 0 
[2] = [Z,11(0)"*] = |-pK 0 ‘ieee 
0 0 0 
Eqn. (2-26) indicates that 
Lge? = opal ee 
Bie ett i,j ¥ 1,2 
i =j : 


(2=26) 


(2-2/2) 


(2277) 


The result given in Eqn. (2-27) is consistent 


Withethiate Guven ineuon.9(2-0)ie fOnusamuUlm_Orlueeatch, 


(2) Layered earth: 


Referring the model to (II-D) and letting the 


primary H-field at the surface be denoted by 


20 


H{0) = MA. 


the total H-field below the surface can then be written 


as 


Hy (r) 


I 
> 
0) 
+ 
ey) 
0) 


Le (4) = Are + Be P 
y y 


Comparing the above expressions with Eqn. (2-16), it is 


clearly seen that -[¢(r)] has the form 


Wee YA le Gerd 
iL x i 
e Tt ehae (= 0 0 
x 
“KZ B Ki2 
Moan) = 0 e + 5X e 0 (2-28) 
M 
0 0 0 
where Ky = Vj wu/ py and Py is the resistivity of the 


first layer, and Boe BY represent the effective reflec- 
tion coefficients from all discontinuities. Since the 
parameters for each layer are constant within the 


layer 


>| 


B 
SS ee 
A A 

x ay 


By following the procedure in the last paragraph, one 


can obtain 
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= 0 0 
1+ 
-1 A 
[>(0) ] = 
ak 
0 
0 B 
Lit x 
e B 
0 0 


and the surface impedance tensor 


0 PK) (ay) 0 
: Seo AB 
[2] = [2,119(0)17*= |-p) x, AGE 0 0 
0 0 0 
(239) 


which is consistent with the result given in Eqn. (2-13a) 


one (2—l0). 
ae tie A - B ey 
ceri en ool ole aeceo (2-30) 
(3) Vertical fault - two dimensional model: 


A vertical fault is shown schematically in Fig. 
(2-3), the x-axis is along the interface of the fault, 


or the strike direction. Py Csi) ee Seethe resistivity 


Crea LoeCm@ GOO LOne a. 


Let the primary H-field at the surface be 
> 
gM AU, ce ERG Oe 


The total H-field can then be written 


-K.Z 
ak 


A, (e om RL V¥r2)) 


Hy (r) 


T 
a 
Oo 


Hh a ga SD 


H, (x) = a 


where R's represent the normalized reflection coeffi- 
cients and depend upon y and z in this case, The 


associated functions are 


-K.Z 
e wenn R 0 
XX 
“Kz 
= 0 + R 2-31 
One) al e vy ( ) 
0 Rey 
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ef T+, 0 
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0 0 
ay Roy 0 
0 0. (K R + _R 
it ee Zee S/n 2, 
[2] my CREE me ae . 
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where sf om = ac [Ruy (2) 


and the surface impedance tensor 


0 Ps 2k y* y®z ) 0 
= Sree aR 
yy 
OAC uate Ny) 
[Zz] = Spy 0 0 (2-32) 
xx 
0 0 0 
Pi {Ki 7 2Byy t yRzy? 
op oe re ae Coney 
yy 
i AWG, Se 
ee ee Zaxk i: 
ees aeRO), e289) 
XX 
Zi, = 0 i,j # 1,2 (2533c) 
il. = 3! e 


The result 2 * shows the anisotropic behavior 


ipa Seal 
of the impedance tensor elements in the two dimensional 


case. 


II-G. Reduction of a Three-Dimensional Problem to a 


Formal Two-Dimensional Formulation 


The following analysis only applies for the case 
where the principle directions of resistivity are 


mutually orthogonal and thus a cartesian coordinate 
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system is applicable. The 9 impedance coefficients 
shown in Eqn. (2-24), are reduced to 4, generally 


non-independent, coefficients. Amongst all the com- 


ponents of the fields, only ix and a are independent, 


while the others, Hoe Hye Hoe Ee EY are linearly 


related to the two independent components, The math- 


SMa eucamepLoguCanmbe, Carried OuLsDY rewriting 


Eqn. (2-16) 
Hye ty ony (2-34a) 
Hy = doiH + bool, (2-34b) 
H (2-34c) 


a a 378 Dx + 32H ; 


It is obvious that Hy, Hy and Ho are not independent; 
and each one of the three components are dependent on 
the other two... To put He in terms of Hy and ake £Or 
instance, H - and H can first be obtained in terms 
(ope Hy and Hy from Eqn.” (2-34a,b), and then substituted 


in -thegkqny (2-34c) 
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oH + 0H, (2-35) 
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An alternative proof can be achieved by consi- 
dering a model in which the “earth resistivity p(r) 
approaches some constant value in the region near the 
surface. At the surface, J, (0) = 0, and thus E, (0) = 0. 
Substituting E, (0) = 0 into Eqn. (2-24c), an equivalent 


relationship to Eqn. (2-35) is obtained, 
H_ (0) = a. HY (0) + py A (2-36) 


Substituting Ho Eromencn. (2-95), On 2—360)=1ntoy (2=24)), 
=e : 
SMM ea lel Bx val mel OLY, 


E = ony 


i] 
eae DCs 


where 


For convenience, the unprimed Ae with i,j = x,y 
Will be used in the following téxt for surface impedance 


for two-dimensional problem, unless otherwise indicated. 
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Hee = eH net Hi (Qeo7a) 
Ea OH recor (2-375) 


The vertical component of magnetic field Hy is a 
secondary field, and is related to the total H-field 


components by Eqn. (2-35) or (2-36). At the surface 


H(O) = o.Hu(0) + ao (0) . 


The relation between HH, and the E-field at the surface 


can be written 


H, (0) = YE, (0) au yBy 10) (2-33) 


which defines the surface admittances [Y,]. Eqns. (2-35), 


(2-37), and (2-38) give the following relations 


Ouee= = Geek ar NAL YE (2—39a) 


QS Y 4 steams vier ° (2-39b) 


As a function of frequency, the impedance and 
admittance are related to the field components ata 
given point at the surface for any polarization of the 
source and for any arbitrary earth conductivity distri- 
bution, and hence they contain diagnostic information 


for geological interpretation. 
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Fig. (2-3). . Vertical fault. 


Fig. (2-4). 


Rotated coordinates. 
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II-H. Orientation Behavior 


Surface impedance and admittance elements aige 
Yay and ot which are functions of the earth parameters, 
contain the information for determination of geological 
substructure including their directional behavior. 
Taking the 2-dimensional problem as an example; as 
shown in Eqn. (2-33), if the axes of the coordinate 
system coincide with the principal directions, i.e. 
the strike direction and the perpendicular, ae and 
At vanish. Physically, in either of the principal 
directions, a component of the E-field is correlated 
only with its orthogonal component of the H-field; 
CHUSmCNem I LNGUCtTOn. parameters Zx and Sony are zero. 
However, if the axes of the coordinates lie in any 
other directions, where the effective resistivity is 
no longer isotropic, each component of the E-field is 
induced by a Combination of all,.components of the H- 
field, and the magnitude of each a55 is a measure of 
the induction effect between the corresponding E and 
H components. 

In order to examine the directional behavior of 
gtsy onemcannrotate the: fcleldsul pdndeleal bOmecilweatLol— 
trary direction x-y, i.e. the measuring direction, to 


a newadirection x"=y'¢as shownein (Rigs) (2-4) 7 where 
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cos@ EB) i+ ¢sine’E 
x ay 


ea] 
II 


-sinOd E + cos@é E : 
x yy 


iniiusedetines the 72Otati1on Matrix which 1s orthogonal 


cosé@ es 


[R] = (2-40) 
—-Ssindg cos6@ 
Orhtogonality of the matrix enables one to write 
[ES] =VER] [5] 
PH Lea Rett] 
[Eee] eZ." pe [HS] 
and the rotated impedance matrix is obtained 
Benne 1 
Pie RIM PsIelR] : (2-41) 


Thus, the rotated surface impedance elements ae are 


related to those in the measuring direction 


Oy al - ee P 

a (owe 2) cos20 + (2, +2.) sin2e] (2-42a) 
1 oli a 2 a z 

TR pg gh acara ead C= Z.) cos26 (2,42...) 8in26] (2-42b) 
’ ik : 
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Sewage lie » = 
Z me ya elit WA ann eee re (2, 


yx xX 
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where 2; (8) aremLunectrons tof ‘the rotation angle 6, 


and also have the following relations 


Z sy (9) SAD ARETE, (2-43) 
Zo (8) = ~Z, (8+90°) (2-44) 
ae EL eee aN (2-45) 
ANG Re Rea Ale Bis (2-46) 


where ay and Z. are new constants, As shown by Sims 


WUJZo9jme che loci of Z2..(0) in the complex plane are in 


q 
4.3 
general elliptical as shown in Fig. (2-5). -This 
graphical representation is very helpful in visualizing 
the impedance behavior as a function of rotation and in 
the determination of dimensionality. 

For convenience all ellipses of ae Der ie 
(2-5) are shown to the same scale and correspond to 


the same angle oe between the measuring axis and the 


Major axis of the resistivity structure. 
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The angle of the major axis can be found by maximizing 


either 


7 v v 
JZ py (8-2) | or [Zjey (0) 42). (8) |. 


d Pa ] xe 
AB | Zacy (8) + 25 (8) | =—07 3 (2-47a) 


The 4 possible solutions of Eqn. (2-47a), corresponding 


to the major and minor axis intercepts, are 


2 Re[(Z_-z__) * ( Zeca 
(2-47b) 


This result is ambiguous since either of the principle 


Al 


directions of resistivity is found. However, this ambigui- 


ty can be removed by apparent resistivity information. 


The rotated admittance elements can also be 


obtained 
i] i 
Y,, (9) = Y, cos@ + yy sinég (2-48a) 
q 1 
Y = -Y cosé@ + Y_cosé6 (2-48b) 
y 69) cs g . 


and are related to 
v [a ’ 5 
Y,, (8) = oe ae 


Thee LOCUS son ¥; (8) is also an ellipse as shown in 
Fig. (2-5). The ellipse is centered on the origin and 
oriented at the angle Oe The rotation period for 
this. loctsyis S60°%1n 6, ansteadvofels0tses for Zia 


THe MaLOreanduligoreaxes are 


’ ‘ i] pal 
major axis 2¥, (6...) (2-49a) 


° 2-49b 
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Fig. (2-5). AN Y; (8) and a; (@) TOC iN 


complex plane. (0=0 corresponding to the measuring 


axis not shown). 
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S15) 


Using the procedure as used in dealing with Dawe 


the angle a Seu. 


e 
PeRe WY) 
mal -1 x 
oe aay tan F z . 5 . (2-50) 


The rotated Ol behaves in a similar fashion to 


' ' ' ' Oar; | 
the Yio except that Then yy and Oy Say 


II-I. Special Cases and the Indications of Dimensionality 


G1) One-dimensional case: 


When the earth conductivity is uniform or is a 
function of depth only, the surface impedance is in- 
variant with rotation angle 6. Eqns. (2-27), (2-30), and 


(2-42) give these results 
Z =a a7, = © (2 ch) 


' 5 pa 
Bras = TN 2 (2-51b) 


where Zy is a complex scalar for all 6. 

The rotation loci for Bias as shown in Fig. (2-6), 
reduce” ‘tompoints in the complex planemcentered zat +25 
(] r] c . i] i] 

Z batel 74s 
for 2 xy and Sox! and at the origin for eae vy 
In the one-dimensional case, no vertical magnetic 


field exists (HL = 0), so that ge = Ys =)0. The rota= 


tion loeL, tox a, and OL are also point ellipses 
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cenveredvon, thesorigin) as Shown in, Fig. (2-6). 


G2) Two-dimensional case: 


As illustrated in Eqns. (2-33) for a two-dimen- 
sional model, Zoe and ae vanish along the principal 
abdei, bain i) = one Oi Os Eqn. (2-46) also shows 


that. the equation, 2% .+2' 


oe Zoy = 2Z,, 1S invariant with 


rotation angle 6, so that, for a two-dimensional model 


Zax ayy = 22. = 0, and thus, 
i] peace i] 
Z ex = ayy Lon-all 6 (2-52a) 
: = : = ao ° = 
Z ex Vey Omron d 64" 6,+90 : (2-52b) 


The loci for Base and are are straight line ellipses 
Geneercdwonsthesovigin, The loci for Zi) and ZN are 
alsopstraighit line ellipses centered at +2), as shown 
say SPeMei, OS 

A vertical magnetic field He exists when the 
primary magnetic field has a component perpendicular 
to the strike. Elements of ve and Oe repeat for each 
180° increment in 6. Suppose y'-axis is perpendicular 


tomthemstrukem(renercings tO Plgs (2-5 Mande 24) je, 


v 
Y, and qa! are 
7” fy 


Van | a. — @) 6 


° ° — 
ea hy (Gyr 0 cme 0o) (2-53a) 


y' =q ae () ts) 


¢ re (ooa 1 B0es) (2-53b) 
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Bugs 2-0 )es Bg hed 4 Xj (8) and as (8) elehh Rene 


one-dimensionality. 


5 


Yj. (8) 


Y¥' (0+90 
y' ) 


Fig. (2-7). rR Yj (8) and 


two-dimensionality. 


ay (8) 


ay (8) 


a’ (6+90 
o ) 
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s) co eile (2—53¢c) 


ze a Ooh 


The loci of ¥, (9) and a; (6) are straight line ellipses 
cencered Ongtnesorigan as shown invE1g.. (2-7) . 
The direction corresponding to the angle oe 


determined by Ze can be one of the two principle 


j 
directions. The angle Oe! however, determined by 
a (otad eee by, a5) corresponds to the direction 

parallel (or perpendicular) to the strike for a fault. 


(3) Three-dimensional case: 


Arbitrary conductivity structures could produce 
the elements of impedance and admittance, ais and Yar 
having arbitrary values, and the rotation loci may 
in general be ellipses as given in Fig. (2-5). How- 
ever, if there exists a plane of symmetry through the 
measuring point, the behavior of aa 5 and Y; with 
respect to the rotation angle 9 appears as those of 
a two-dimensional model. These special cases of three- 
dimensionality can be classified into the equivalent 
two-dimensional problems; some examples are shown in 
Fig. (2-3) thessimplitication 3n) thesescasessresults 
from the loss of information caused by making measure- 
ments at a point of symmetry. In fact the generalized 
two-dimensional impedance Zi and admittance Y, can 
lead to unique interpretations only for two or one- 


dimensional models. However the more general result 
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Fig. (2-8). Three-dimensional models with a plane 


of symmetry, through the measuring point @. 
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can give information on the dimensionality of the 


Structure in terms of the dimensionality indicators 


‘a Zt Ze t% ‘B 
a (skew) = ge = poy = uEr (Swift 1967) 
KY ey x Vay x 1 
(2-54a) 
J q 
=, 
6,2 (Word et al. 1970) (2-54b) 
zy ee 9 
yo yx|'0° 
nee 7s 
oo 2 (Word et al. 1970) (2-54c) 
Van Ey 
XY yx os 
and 
y! a 
Boe=— ee ents (Word et al. 1970) (2-54d) 
' zh a 
vy ) Y Goh 
ze 


Skew ao is invariant with rotation angle 6, and 
equal to the ratio between the two displacements from 
the origins to the centers of the corresponding 
ellipses as shown in Fig. (2-5). It is obvious that 
the condition a = 0 can indicate either one or two- 
dimensionality. Word et al. (1970) propose that a=0 
is a necessary but not a sufficient condition for two- 
dimensionality. The sufficient condition must include 
the condition Bo = 0. Bo is the ratio of minor to 
major axes of the Zi 5 loci as shown in Fig. (2-5). 


For two-dimensional cases or equivalent, the rotation 


rt 
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: Te 
é) y 70h ae 4 
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loci are line ellipses, so that BY is zero. However, 
for the one-dimensional cases, By = (0/0) becomes an 
undetermined quantity. 
For one or two-dimensional cases, Bwe = 0. shi 
has the advantage as an indicator by avoiding Oe ee 0/0 
for one-dimensional cases, in addition, it also has 
potential for judging the degree of three-dimensionality. 
B = Bo, = 0 for two-dimensionality, but is not 
necessarily non-zero for three-dimensionality. 


The behavior of the dimensional indicators are 


summarized as following: 


One dimensionality a = 0 
So hee 0/0 
Boo * . 
Bre z Boh cece 
Two dimensionality a = 0 
and special ee ame 0 
case of symmetry cae = 0 
Boe st Boh att 
Three dimensionality a * 0 
sa oe 
Soe, te 
(may ) 
Bre * Ban = Our 


yh 


6 oe ata 


a ow 
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While the use of the dimensionality indicators 
provides some information regarding the nature of the 
structure, quantitative interpretation is restricted 
to one and two-dimensional models. In practice it is 
frequently possible to interpret the more usual three- 
dimensional structure in terms of these two and one- 
dimensional models to some reasonable degree of 
approximation. In this case, it may be desirable to 


define the following parameters: 


q 
ehepeay' Je 
8 2 (2-55a) 
Pay * Pyx 
q i] 
Py ~ 
Bo = iF (anisotropy factor) (2-55b) 
PM 


where P44 are the apparent resistivities corresponding 
' i] i] : 
y um and 
xe) ay to)? and py and op, correspond to the maximum a 


i] 
minimum values between oe and Pyx' 


3-dimensionality By * 0, measuring the degree of 3- 


dimensionality. 
B, = 0-1 
2-dimensionality pe 0 (ome = Peng = 0) 
Soy te WE We Gea! 
l1-dimensionality By = 0 
Boks Clg OL a! 0.) 


CHAPTER VELL 


INSTRUMENTATION AND DATA ANALYSIS 


IfII-A. Instrumentation 


The instruments used for data gathering have been 
described by Allsopp et al. (1974), andare shown in 
Fig. (3-1). It is designed to measure three compo- 


nents of the magnetic field (HL, H 7H) and two 


y, 
components of the electric field (E/E) of the 
micropulsation signals in the frequency range from 
OR 000i co, LO Ho: A short description of this system 
will follow. 

(1) Magnetic system: 


A time varying magnetic field with amplitude of 


s to nos gammas is detected with a 


ehemordensot 10m 
specially designed induction type sensor (Burke and 
Rankin 1975) and produces. uV level output in the 
frequency range of interest. Preamplification stage 
consists of a RATEK RA-1l parametric amplifier and two 
FET input operational amplifiers, with a total gain of 
85 dB. The output is fed to a post amplifier with swit- 
chable voltage gains at -20, -10, 0, 10, 20 dB, and 
then through a Butterworth low pass filter to eliminate 


the aliasing effect. The overall amplitude and phase 


responses are shown in Fig. (3-2). 
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Fig.(3-2). Amplitude and phase responses of the 
Magnetic and electric systems. (1) for low mode and 


(2) tor high) modesor recording. 
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(2) Electric system: 


The electric field variations are detected as a 
voltage difference between copper electrodes in the 
ground. A chopper stabilized operational amplifier 
is used for preamplification. The post amplifier and 
aliasing filter are identical to that of the magnetic 
system as described above. The overall amplitude and 
phase responses are shown in Fig. (3-2). 

(3) Digital wecording: 

The recording system consists of the multiplexer, 
A-D converter, crystal clock and a seven track tape 
recorder, with a dynamic range of +10 V. The sampling 
rates were chosen to be 1.25 and 40 samples/second/ 
channel in low and high mode respectively. A 2700-feet 
PApeCHWECI MoD Om bp. Call last up to 20 days or 16 hours 


depending on the mode of operation. 


III-B. Data Analysis 


The first step in the interpretation of the data 
is to obtain power spectral estimates of the measured 


quantities, Hy H 7H E E 


y Xsan Vow 


(1) Power spectrum computation: 


Let f(t) represents the time sequence of a re- 


corded signal. The procedure is as follows. 


* 3 ae 2 ¥ 
Lvehen Ona Ro esos 
} wo =p oe ‘Gvia 
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Bm Cuncation 

Let N represent the number of samples per data 
set; To = NAt is the sample time length (sec) per 
data set, where At is the sampling interval (sec). 
A data set with length of N samples can be obtained 


by truncating the possible infinite time sequence f(t) 


£. (t) = £(t) W(t) (SE) 
where 
Ts roiqe MOLE aS twas T)) 
Wy ft) ss ° 
0, (t = otherwise) 


The effect of this convolution is removed in the fre- 
quency domain. 
b) Removing mean and trend 

Since only the AC component of the signal is of 
interest, it is customary to eliminate the DC compo- 
nent at this stage by removing the mean. Removing 
the trend is done to remove the effect of contributions 
at periods longer than the fundamental and is accom- 
plished by removing the slope corresponding to the 
difference between the average of the first and last 
Lil ras pOLnes. 
c) Tapering 

A large discontinuity may occur between the value 


of the function at the beginning and the end of the 
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data set; this discontinuity produces undesirable high 
frequency harmonics in the power spectrum. It is 
customary to taper the data set making the first and 
the last points approach the mean value which is now 
zero; and it is done by using a cosine bell window 
with one-tenth of the samples tapered at each end of 
the data set. 
da) Fourier transform 

The Fast Fourier Transform method is used to 
transform the time sequence F(t), from the time 
domain into the frequency domain. 


2 f° 
Fo.) = ge |. fy (te 
* 0 


ST 
— * = 
dt F (yw) Wy y) (3=2) 
where F(w) and W, (w) are the Fourier transfer function 

Om f(t) and Wy ft) respectively, and w = 27f. From N 


data points, N/2 complex coefficients F (wy) are 


obtained with 


Wy =Nw, ¢ where W) = 2n/T MAT COMET ean yp otetes ep N/a 


The corresponding frequency range is 


(£,=1/T.) < fn < (£=1/2At) . 


The highest frequency fy is called the Nyquist or fold- 
ing frequency. Components of higher frequencies than 
fy are analog filtered out of the recorded signals, to 


avoid aliasing effects. 
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e) Hamming window 

The Fourier coefficients Fw) obtained from a 
truncated data set are different from the actual ones 
F(w) due to the truncating effect. The correction 
can be done by multiplication with W, (w) in the fre- 
quency domain. However, the Hamming or a similar 
window can be used as an excellent approximation for 
this correction. The Fourier transfer of the Hamming 


window has the form 


W, (Wy) =0.54 Wi (w,) + 0.23(W (way) + Wo (wy). (3-3) 


Applying this window by convolution in the frequency 
domain, the Fourier coefficients of the data signal 


are obtained as 


F(w,)= 0.54 Fo (w,)+0.23[F (u,_3) + Fo(wi yy]. (3-4) 


n+l 


£) System response correction 


The overall system response must be removed to 
obtain the true Fourier coefficients of the E and H- 


field quantities. 


g) Power spectrum computation 


Auto and cross-power spectra are computed as 


— * -_ 
Prz(w,) = Fy (w,) -F; (w,) (3-5) 


where * indicates the complex conjugate and I,J 
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(=1,2,3,4,5) are integer indexes referring to the data 
channels. The number of samples per data set used is 
usually large (say, 4096), so that there are a large 
imber (N/2) of harmonics in the power spectrum. At 
this point, smoothing is carried out to obtain 
average power spectral estimates. This reduces the 
computing time in later stages of the analysis, and 
eliminates noise effect (see Chapter IV). Ten fre- 
quencies are obtained in each frequency decade; the 
ratio of center frequencies of successive bands is 

f 


k+1 
fy 


= 

10 > 

Thus, the bandwidth averaging power spectral estimates 
at a representative frequency Wye is 


aL 


Tea 1 A) (3-6) 
b (band) 


a es Te on 


<Pryz (wy, 


where M,, is the number of harmonics involved in the 
bandwidth averaging. 


If more than one data set is used, then the over- 


all averaging spectral estimate is computed as 


M 

s 

Ss m Ms 
s m=l 


where m = 1,2,...M, represents a given data set and M, 


is the number of data sets contributing to this average. 
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(2) Polarization: 


The polarization characteristics of the MT-field 
have been studied, using power spectral estimates, by 
Fowler et al. (1967), Paulson (1968), and Rankin and 
Kurtz (1970). The three important polarization para- 
meters are: Degree of Polarization, Polarization Angle, 
and Bllipticity. 

a) Degree of Polarization (DP) 

For quasi-monochromatic EM-signals, the power 
spectral estimates may consist of polarized and un- 
polarized portions, and can be represented by a matrix 


of the following form 


= + (3-8) 


where x and y refer to the components of either E or H. 
Thestirstematrixyon the right hand side corresponds to 
the completely polarized contribution, i.e. that part 
with unit coherency (AC= B*B); while the second matrix 
corresponds to the Deen erantecontemormens The degree 
of polarization is defined as the ratio of the polarized 


intensity to the total intensity of the signals. 


youaae too ‘of fn wu 


3 
ete ar ae ai 


or 


7 
ee on 
7 —_ = |S > 
or ae. = NCL ot Ee cub 
ay 7 1 ¥ 


Sp! 


fy A+cC r 4|\a| 2 
re Sx cities: a i ne (<P) Sort < 2 ae 
vy, Bese Pyy”?) 
where J = < > - i 
ee seyy” 5 Pay? <P yx?! is the determinant 
of the generalized power matrix. For a completely 


polarized signal, une coherency between x and y com- 
PONenL Sei Seunityecomthat, Ja=0 and DP = 1. For a ton- 
pletely unpolarized signal, where the coherency is zero, 
DP = QO. 
b) Polarization Angle (0,,) 

The polarization angle 85 is defined as the major 
axis azimuth angle of the polarization ellipse of the 


polarized field 


alee = ° Set 
85 ( ) 


Coe Cae) 


Ellipticity is defined as the ratio of minor to 


Major axis of the polarization ellipse 


H 2 Im<P. > 
¢ = tan|% sin 7 


a % 
Ie <Pyy?) + 4<P Ly ><Pyy?! 


(3= 181) 
€ values ranges from 0 for a linearly polarized field, 
COml~mLOm aecircularlyepolarizedut beld.mmA, posi ti Vvesor 
negative value of e« represents the sense of polariza- 
tion of righthand or lefthand as measured when looking 


into the propagating wave. 


(3) Impedance computation: 


In order to determine the surface impedance 
matrix elements iar using Eqns. (2-37), two independent 
data determinations required. However, determination 
can be done by using power spectral estimates in the 
frequency domain. Consider the frequency domain 


expression defined the surface impedance tensor, as 


ih oh AL le Se, 7 Alem (3-12) 


[ (2-37) ] 
E 


CB Se Yes 
v4 yx x ey ay 


Let each of Eqns. (3-12) be multiplied in turn by the 
complex conjugate of each of the 5 field quantities 


(HLH rH EB). If the source field is randomly or 


yi 
partly randomly polarized with respect to the frequency, 
there is some degree of statistical independence between 
the cross-powers of the various field components. 
Consequently, the power spectral estimates <P> 5 (w) > 


as determined in Eqn. (3-6) will produce 5 linearly 


independent equations fLor each of Eqns.) (3-12) of {the 


form 

as <EXE.? = Zs SESH? ate 44 SESH” 
Zs eee = 25 “EH? aa 24 Ee 
ve <HYE,? = Zs SHY? of 25 HSH? 


a2 


a] 5) 


4. <Hi hee mee cH AHs > bh O7. <oXHS 
Vel vee Nes LV Vey: 
. _<H*E. a ¢ * ’ * = 
5 H*E| > Z 5 SHEH, > + Bessthsti (3-13) 


where i = x,y. 

Using Equations 3 and 4 of (3-13) is equivalent to the 
least square error fitting method (Sims et al. 1971). 
Eqns. (3-13) provide 10 pairs of simultaneous equations 
for solution of each impedance element. However, the 
nee (7 ane ceo) alia Lhose involving (5) are not 
Suitable since under some circumstances the solutions 
for impedance functions become indeterminate. There 
ibes4arcemainiungqupairs. (1/2), (1,3), (2,4), (3,4) which 
are well behaved as long as the source field is at 
least partially unpolarized. Thus, up to 4 indepen- 
dently computed values of each of the impedance elements 
Z.. can be obtained at each frequency Wy It is useful 


1) 
to average 


L 
tL & 
Z.« (w,) SS d PS AG (3-14) 


ug Q 


where 2 = 1,2,3,4 represents a given solution and L is 
the number of acceptable solutions among the possible 4. 
In order for a solution to exist by the method 
described, it is necessary that the smoothing be 
carried out i.e. at least over 2 frequency estimates, 


or to use 2 independent data sets. On the other hand, 
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in order to obtain satisfactory value of Zeay the 


J 
bandwidth used for smoothing must be limited to a 
range for which Ane does not vary significantly with 


respect to frequency. 


The admittance elements defined by 


Lee Ore ele El (3-15) 
[ (2-38) ] 


and the coefficients a and Oty in the relation 


Heea= 10 Bn chao) (3-16) 
Zz opr uy 


[ (2-35) ] 
are computed in a similar manner and are subject to the 
same conditions as the impedance elements ai5 just dis- 


cussed, and the result is 


L 
a Q 
Vee ae Ye (Gos) (3-17) 
1 L 2=1 cee < 
L 
ae v Q 
Cha Eze yh oie AUR DP BA (3-18) 
1 L Q=1 a a 


The apparent resistivity for an impedance element 


ais is computed as 
= 22 Zea : 3-19 
P44 (Wy) = F [255 (4) | (in ohm-meter) ( ) 


where the unit of a is (millivolt/kilometer) /gamma. 


aD 


(4) Criteria of data pe lvability : 


Experimental data is inevitably recorded in the 
presence of noise of one kind or another, and the 
establishment of reliable criteria for judging the 
quality of a computed result is important. For the 
case at hand, noise is defined as any portion of the 
measured EH-field quantities not conforming to the 
impedance and admittance relationships defined in 
Eqns. (3-12), 43-15), and:.(3-16).. Both phasor 
coherency and predicted coherency, have been proposed 
and used as a measure of data quality. 

a) Phasor Coherency (CP) (Word et al. 1970) 


When noise is present in the measured EH-fields, 


the independently computed solutions of Liye ee and 


ay are not completely consistent. The degree of scatter 
between the independent solutions can be used as an 
effective measure of the reliability of the results. 


The phasor coherency, representing a measure of the 


degree of data quality, is defined as 


CP = 


: 


where R 
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and Rp represents the %-th solution for 7h cer or a, 


and L the number of acceptable solutions. The phasor 
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coherency has the range of 


Ly Seis 2k, y 


bm PredivcteamCoherency @(PCH) (Swift 1967) 


The predicted coherency is defined as the coherency 


between the measured field component and the correspond- 


ing component predicted from the average value of Lise 


ae or ay. The PCH for the impedance is given as 


[Zee Beno yt eS E XH | 
PCH (E*EP)= pe xX 5 X Xx 
x X r, 
< > < >+ * <p* > 
t<Pi? (2,5 Prax alae Puy Ren ay Ay Mh 


DEC ESTE tee eH > | 
PCH (E*EY)= x Z 
{<P..> (12, | 


26) 
<Bap tl2yy 


SPE >+2Re(Z* Z <H*¥H > 
py T2RAZS 2 HEH 7) 


and for the admittance 


(ia SE RANE oae rn athe > | 
Pou (ig) = =. : : 
<P_>+ STOR Pe 
pee tire <P, ty, Paap t 2Re(Yy Y <EXE,, )] 


where the superscript p represents the predicted quan- 


tity. The PCH has the range of 


oO 
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PCHS<ee] aes 
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The phasor and predicted coherency represent 
a measure of the degree to which all influences in 
the measured EH-fields are related by the same tensor. 
Gees leands ? CHe—eleiteall solutions of a given Re 
coincides in both magnitude and phase. But the con- 
verse is not necessary, since noise can exist and be 
coherent between channels. As part of the process of 
analysis, the ordinary coherencies (EH) and (EH) 
are computed for the raw data. If these ordinary 
coherencies are low the data is rejected from the 
subsequent processing stream. 


The problem of noise contamination will be dis- 


cussed in the following chapter. 


CHAPTER IV 


NOISE DISCUSSION 


In the statistical sense, the primary magnetic 
Signals appearing in the magnetotelluric method can 
be considered as noise. However in arbitrary measuring 
directions the orthogonal components of the magnetic 
fields, consisting of primary plus secondary contribu- 
tion, show a considerable degree of coherence. The 
components of the electric field at the surface of the 
earth are substantially secondary and are thus highly 
coherent with both the orthogonal magnetic component 
and between themselves. In the subsequent discussion 
those components of the fields which are coherent in 
the sense just described will be considered to be signal 
whereas all other contributions, whether coherent or not, 
will be considered as noise. 

As an example of coherent noise, the 60 Hz power 
line contribution is frequently present and strenuous 
efforts must be made to remove it. Other examples are 
power line spikes and steps, which can affect all 
channels of information. This noise is usually elimina- 
ted by discarding sections of record during the pre- 
editing processes. 

Random noise cannot be eliminated by any process 


except in a statistical sense, but since it is incoherent 
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its effect can be taken into account and noisy data 


can be rejected according to a predetermined criterion 


for coherency. 
IV-A. General Discussion 
(1) Covariance and Coherence 


The covariance functions for time sequences can 


be defined by 


a,-(T) = a ) X(t) X,(t+t) auto-covariance (4-1la) 
ik Ni eH1 ib Hy 
1 N=-T 
Crzy(D RF a X, (t) X, (ttt) cross-covariance (4-l1b) 


where Xr and Xs are time sequences with zero mean values, 
aiding OG AZgi. oa eeLepresents!-the time lag. 

For signals which are quasi-periodic and coherent, 
a(t) and c(t) are also quasi=periodic functions of T. 
For random noise which is non-coherent, a(t) and c(T) 
are also random, and when N is large a(t,tX0) and c(T) 
becomes very small and can be considered negligible when 
compared with a(0), where a(0) is the zero lag auto- 


covariance which is also the energy density of the time 


sequence. 


N 
a (0) == ) Ix(t) |? = 9° (4-2) 
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Fig. (4-1). Relative amplitudes of auto-covariance 
a(t) and cross-covariance c(t) for two coherent 
Sinusoidal sequences with a phase 9 (upper), and 


for two random noise sequences (lower). 
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This also turns out to be the variance estimate of X(t) 
which, if X(t) is stationary, approaches a constant 
value when N is large. 
Fig. (4-1) shows the relative amplitudes of a(t) 
and c(t) for two coherent sinuisoidal sequences with a 
phase difference 6, and for two random noise sequences. 
Now consider time sequences which consist of the 


sum of signal S(t) and random noise n(t) 


X,(t) = S(t) + n,(t) (4—3a) 


X5(t) = S, (t) + nN, (t) 5 (4—3b) 


The auto- and cross-covariance functions are 


1 oe 
a_(T) == Syl) 3S (mi aes) 
I Nettie 1 I 
1 {Net N=-t 
== L ae tt ee ee ae) 
Se AG | 
+ S(t) n- (t+T) + n-(t) S_(t+T)}] (44a) 
22 eee eereeear) Ot 
N-t 
c_(t)== J x. (t)X_(t+T) 
85) ee eae a 
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= DE Se es ae AGS US) 
ay N-T | 

fe S_(t)n_(ttT) + N_(t)S_(t+T) (4-4b) 
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Inewiinstacermeinebothwa (mtandec(s) i1sesirgnal 
and the second term is noise; the last two terms are 
the noise-signal interference terms and pertain to 
the noise effect. Except for the second term of the 
zero lag auto-covariance a Oe which is also the 
noise energy density, each of the remaining terms con- 
taining n is a sum of random functions as indicated in 
Eqns. (4-4), and thus decreases as N increases. For 
sufficiently long time sequences, it will be shown 
later that terms in Eqns. (4-4) involving random func- 
tions can be neglected. 


(2) Power Spectrum of Signal and Noise 


a) Power spectrum of random noise 


For a random time sequence n(t), its complex 
Fourier coefficients are defined by 


=Jot 
nw) = - [ mictoe edt (4—5a) 


0 


which for digitized data becomes 


> N -j27mnt/N 
T(r es eee) ert) : (4-5b) 
=], 


The corresponding power spectrum is defined by 
= * = 
P (Wy) n*(wi)n(w) - (4-6) 


n(w) and we are also random in both amplitude and 


2 ie 
te 


7 / _ 7 
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phase. For digitized data one obtains for the average 


noise power 


where 


> N/2 
N dL Se ny 
N 
ee B ie mia i 
N v2 ES) ee on 
N 
== J |ntt)|* (=a(o)) 
t=1 


(4-7) 


is the variance estimate which approaches a constant 


value as N becomes sufficiently large. 


Henceforth we 
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will consider discretely digitized data of finite sample 


length N. 


As mentioned earlier the distribution of the 


amplitudes of a random noise power spectrum is also 


random. 


According to an experimental study of a 


typically random process described in section A-4 of 


this chapter, the maximum value in the spectrum 
Fw.) ~ 4€. It can be considered that the power 
spectrum Sone can be represented by the average 
value € weighted by a factor K* so that 


Pb@n? = K 


The Ke form a set of random positive numbers with 


N 


Z 
n 


<P ee 


(4-8) 
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Now let the weighting factor for the maximum value of 


P fw) be represented by Kee The magnitude of Ke is 
found to be limited; for instance, ug ~ 4 for the case 


M 


given above. Also, as mentioned earlier, ae has a 
constant value, and thus it is clear that for a sta- 
tionary random process the amplitude of the power 
spectrum as given in Eqn. (4-8) is inversely propor- 
tional to N, the number of sample points. 


b) Power spectrum of signal 


For a monofrequency sinusoidal signal S(t) = 


S,Sin(wot), the Fourier coefficient is 
T 
= 2 | 2 
S(w,) = 7 | Ss sin(w tie dt (4-9) 
and the magnitude of the power spectrum at this fre- 


quency is 


P.(w,) = S*(w)S(w,) = so (4-10) 
which is independent of the number of sample points N. 
The same condition holds for each harmonic of a more 
general time sequence. 

It is seen from Eqn. (4-8) and Eqn. (4-10) that 
an obvious distinction exists between signal and random 
noise. For a signal the amplitude of the spectrum is 
independent of N, while for random noise, it is 


inversely proportional to N. This fact can also be 
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understood by considering that for a continuously 
periodic signal the power is distributed to only 
those frequencies that are present in the signal, 
while for random noise, the power is distributed 


over all N/2 harmonics with equal probability. 
Cjeepectral Censities withy mixed signal and noise 


For a mixed signal and noise time sequence as 


given in Eqn. (4-3) 


SC) pero GE) nt) (4-lla) 


the Fourier coefficients are 


X(w) = S(w)+n (w) (4-11b) 
where . 
2 -iw,t 
S(w) = Z| S(t)e at 
0 
and 
Je , 
2 a 
n(w,) = = | ni(t)e oer 64 
0 


The auto-power spectrum is 
= = 4-12 
P (uw) = X* (w) X (w,) Pera P Wy) (4172) 


; caren 
where the signal is Pow) =S (w)S (wy). and the noise 
contribution is PL (wy) = 0% (wy) n (wy) +8* (ny) n (wy) t 

n* (w)S(w,). Using the result deduced in the last 


section, 


Soe Sore 
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(4-13a) 
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iL Pig Pi ik 
7 40 Kn + — 4s ,0K,,cos (y,) (4-13b) 
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where 1, is the phase angle of S*(w)n(w) and is random, 


ranging from 0 toamT. 
The noise-signal ratio of auto-power (n/S) . for 


a given harmonic is defined as 


P_ (Ww) 2 
a Pepe Sha) 1 40 2 1 4¢ 
AS) poe == K” + — — K.cos(y..) (4-14) 
a P(w,) N s2 n JN So n n 

and since Kh Ss Ku and cos (y,) can be either positive or 
negative 

ay aC ae 

(5 Ky Ky) S (0/8), S$ (= Kyt Ky) (4-15) 
Ns VN Sys Nso YN Sy 


The general cross-power spectrum between X, (t) 


and X5(t) for a given frequency wis 


- x* bs 2 
Prg () = X7 (W) 5 (W) = Pays (4) + Pity (4) (4-16a) 
where the signal is 
* = 
Porgy () = S7(w)S5() (4-16b) 


and the noise contribution is 


Pg (@) = nz (w) ng (w) + Sp (w) ng (w) + NEw) Ss (w) (4-16c) 
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which can be written as 


Bie) 2s_0 39 
Kn a a ee e +iK ee ) 


ly 2 
J (Koy nJ 


(w) =W 40 oe 


Pig I 


where 90's are the phase angles for the corresponding 
terms in Eqns. (4-16b,c). 
The noise-signal ratio of cross-power for a given 


Harmonie 1S deminea as 
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It is seen from Eqns. (4-14) and (4-17) that the 
noise-signal ratio in the power spectrum depends upon 
N, the number of sample points used for analysis. In 
some cases, the power estimate resulting from an 


averaging is desired; these will be discussed in the 


next section. 


(3) Power Estimate and Averaging Effect 


In some cases such as MT prospecting, the power 
spectrum which is used for further analysis consists 
of smoothed estimates, averaged over a frequency band 
(Mb), or averaged at a given frequency over a number 
of data sets (Ms), or both. The resulting auto-power 


is estimated by 
<> egcPeo te <DD>, (4-18a) 
S n 


From equations (4-13) 


* Ms 1 Mb 
QB Sy a ee ) —— ) Po ho) 
S Ms ay Mb ea Ski an 
M —— 
il : D z 
== ly Po(w,) = s (4-18b ) 


where P., (w) is the signal power density of the k-th 


Sk 
M 
set, M = Ms*Mb, and ) stands for the double summation. 
n=1 
1 M 
< > SS = 
ES ui el Pen) 
2 M = M 
= a = ) ee pee F ) (K,cos(y,)) . (4-18c) 
n=] YN n=1 
il M 2 M 
= ) K* + 1 as previously shown. Define } (K_cos(y_))= 
M pe me n n 
<K?, Eqn. (4-18c) becomes 
D =" \Kre 
See = oe + 40 ee (4-19) 
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The noise-signal ratio of the auto-power estimate 


is estimated from Eqn. (4-18b) and (4-19), 


<N/S> = = OS TT C (4-20a) 
S 


or, by using the reasonable estimate of akec 


2 40K 2 40K 
gi <ti/S>, (+ (4-20b) 
Sy YN SM Nez WNSM 
The averaged cross-power estimate is 
TT oie wm Lo setae 
where 
M M Ade 
eel: eel Zz sn a 
Pos? =i p Pet en) =a p (sve ) (4—21b) 
n=1 n=l 
2 M 8 
Bd4:0 coer li nn 
Sino ena i vb (TR aa® ) 
- M 379 M 38 
+ 298 = ) Ke BM at ) Kage ie : (4=21c) 
VN geal * n=] 
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Since the signals are coherent, and assuming that e 
are in phase within a given band, Eqn. (4-21b) can be 


rewritten as 


Be) 
Res (aoe 


However, the phases of noise terms are random, so that 


Eqn. (4-21c) can be written 
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SB cgay? SY SE a eT pemee (4-22b) 


where Ss represent the average of the random phases, 
ranging from 0 tol2r7. From Eqns. (4-22a) and (4-22b), 


the noise-signal ratio of the cross-power is 


jn 0n> sone 
<p Ig° 4 De mak > 6 nN <K > e q 
SIJ ae WN S 
S) 
(4-23a) 
407K 40K, 
< (————_ + ———) .. (4-23b) 
o eon YN SM 


Eqns (4-20b) and (4-23b) represent the limits of 
the noise effect in the auto- and cross-power estimates. 
However, both terms of mat ai as given in Eqn. (4-22b) 
and the second term of sea as given in Eqn. (4-19) are 
average values of random functions, and thus become 
negligible for sufficiently large M. In other words, 
the noise contribution can be largely removed from the 
cross-power estimate and partially removed from the 
auto~power estimate. In addition, there are generally 
some mechanisms which can further diminish the noise 


contribution from auto-power estimates through the use 


of a good cross-power estimate (White 1973). 
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(4) Experiment 


For the purpose of investigating the noise effect 


quantitatively, three noise-signal mixed sequences were 


generated. 
X,(t) = S,(t) + nz(t), Tres el oe « (4-24) 
where 
N 3 
SC ee eo esinuir to.) with cus—p 10. 


The three independent noise components nz were taken 
from the random noise sequence provided by the Computing 
Services Library of the University of Alberta. The 
noise sequences are normalized to of =—o\, and Ke ~ 4 

for these sequences. The formulas used for computing 


the noise-signal ratio are: 


For individual harmonics 


Be (0) eae (Oe) 
(n/S) = whe: as (auto-power) (4-25a) 
ib P (w) 
Sil 
P (Ww) - P (w) 
(n/S)_.= a Se! (cross~power) (4-25b) 
153) Pe) 
For power estimates 
< Dave gate D> 
Tt /.o a Se eee oe (auto-power) (4-26a) 
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SL 
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we eX LS SIJ i 7 
<n/S? +5= =p = (cross-power) (4-26b) 
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where 
Py (w) = X}(w) Xz (w) Poy (w) = 8} (w) Sy (w) 
* 

Bg lee XT (w) X5 (w) Pory (w) =8; (w) $5 (w) 
Tae ie 

SPyp? = ig) Py (o) <Poy? = i ) Pgz () 
1 ii he 
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The computed results are given in Table (4-1) 
With various values of M and N. In the Table, the 
first three columns are the noise-Signal ratio of 
the auto-power for the three channels defined in 
Eqns. ( 4-25b) and (4-26b). Data in the line with 
each M=-number are the theoretical values defined in 
Eqns. (4-15), (4-17b), (4-20b), and (4-23b) with the 


same values of M, N and cos (y,) = tlie 


IV-B. Investigation of Noise Effect on Impedance Elements 


(1) General Solution 


The surface impedance elements S55 are computed 
fromthe. 4s paimsali( ln 2 ae (Loire oie aie ee OL 


equations (III-B-3) 
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and 4 solutions for each element are obtained. These 


are represented by 


a 
ip ental (4-31) 
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where £(=1,2,3,4) represents a given pair, and Nis and 
Dy are functions of the power estimates. The 


generalized power estimates are considered to contain 


the signal and the noise components, thus D and N can be 


represented symbolically by 


MED A SOANOR Seb’ S =p eee. 
as 0 ay aes ot Ke) ae 


where Noe and Den are the signal, and ANY and ADS 
iba Ke) ijo ij ij 

are noise components. 
When the noise-signal ratio is small, Eqn. (4-31) 


to Eirst order becomes 
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ys is the major component of the data and 
thus is a good approximation from which to deduce 


uniquely the resistivity structure of the earth. 


ANS Nowe ADS 
Ad, . Se ee is due to the noise and varies 
I aye p*. pe 
Safe. ate) Tkyfe: 


with different £ value, since noise is random. For the 


Case Of W5—" (394) and dT = x; Ni and Pe ae. 


== ns * * 
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ait Pe Pee eee PS <p se cH ese hn oS Le o<H*H >< a 
ay, Ve eT) Va x yn 
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N =<P.><q*p > = <H*E =< H S 
XX Hy > = y x By. 
AN <P ><H*E > Pon oh ><p > =< Ee Sen Le > SoHo S<H'R > 
XX Hy ane. th xe yas Vee x y n ay Voc 
* * * 
= > =_ < > 
Nyy Sage SHE.” a ae HAY 
* & Ka * & 1* 
= <« <a =< >< > -< >< Hector 
AN Pp gE a pe Pax 7 HEY HOH : HH. a 


Xy H 
After a tedious calculation, the first order approxima- 
CioneoL thes Lous. AZ. are written in the following forms, 


where the superscript denotes the pairs of equations 


used from (4-30). 
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Dei = <PH gE yhx re = <EOH > SHH 
Deal 78g ate ei <HOH > <AyHY> 5 
Dieta = <EYH > (<SEQH = <ESH > <EVH >, 
Dias ~ Cha Pas ae i SBE oe Re 


(2ZJelmpedancemoolutions for Noise-Eree Cross-power 


Estimates 


Supposing that noise exists in the auto-power 
estimate only, and not in the cross-power estimate as 
is often assumed, then all terms in Eqns. (4-33) having 
a noise cross-power contribution will vanish. With 
this assumption, and some rearranging to make all the 


D's real and positive Eqns. (4-33) can be rewritten as: 
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If Mis large, the noise component of the auto- 
power estimate as presented in Eqns. (4-34) is biased 
positive as can be seen from Eqn. (4-19) 

Gian s 
402 40s Kn 


<P > = oo (ae (4-19) 
q /N M 


The second term on the right hand side can be positive 
or negative whereas the first term is always positive. 


An average taken from solutions of Eqns. (3-20) or (4-34) 


ie 
ao 1 My i) Fex3) a 
SA,g> = Bi, +e ) (AZ i 5) Se et <2 55? (4-35) 


provides a better estimate. 


As just described, the noise contribution to the 


X 


auto-power estimate causes noise components Adis in 


the impedance elements, and the average <2 5> can 


become more accurate than the original Erte Consequently 
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this improved aie is consistent with the existence of 
auto-power estimates which contain less noise than the 
original ones do. Thus, substituting BBS back to 


Eqn. (4-30) gives the following relations: 
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Solving Eqns. (4-36) for <P> provides improved auto- 


power estimates corresponding to the more accurate Bane 
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A cyclic operation 
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Vie ao’ (4-37) 
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can finally result in solution having the noise substan- 


tially removed. 
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These final solutions are consistent with the conditions 


corresponding to a unit predicted coherency, which are 
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If the cross-power estimates are noise-free, the 
following conditions can be simultaneously achieved: 
a) Predicted coherency is unity. 
b) 4-solutions of zi are identical. 
As a consequence the apparent resistivity curves are 


smooth, 
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Figs. (4-2) show the result of the operation of 
Eqn. (4-37). Data was taken from a typical section of 
MT record. M, at the high frequency end ranges from 
1000 to 3000. Such an averaging may remove the noise 
contribution from the cross-power estimates and leave 
noise only in the auto-power estimates. The 4-solutions 
Sue Zee SrC@tie Geeaveragce-4y .dle ploulea, eLom Lets to 
right in each row. The first row represents the raw 


13 
of cyclic operations indicated on the individual plots. 


Z.. and the second row represents a4 after the number 


The scatter in as 5 appears in the raw data, however, 
after a certain number of the cyclic operations, the 
noise in the auto-power estimates is substantially 


removed, and then the solutions of aT become unique. 


(3) Impedance Solutions as Noise in Both Auto- and 
Cross-power Estimates 


When noise exists in both auto- and cross-power 
estimates, the solutions of ans, are those given by 
Eqns. (4-33), and an average value given by 


O 


SPAR ed iy VAR, ~ ek WA nee 
Le] 1 


J 


-also provides a noise-reduction-effect. Consecutive 


manipulations as performed by Eqns. (4-36) and (4-37) 


results in a set of <Z55>", where Me=— li, Zips 2 eMC 


oe 
I 


the number of cycles performed. None of <Z is 
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expected to be noise-free, and it is even possible that 
the noise content is increased. Whether an average 


given by 


Rai IS potas I A (4-39) 


can provide a noise-reduction-effect depends on the form 
Oteene nOLsewcontcribucion. sHOWeVer, 1t 1s significant 
that even when no noise reduction is effected the 
results can be used to establish a confidence limit, 
which can be used as a criteron for final data selection. 
The confidence limit is defined by 


CL (4-40) 


II 
Be TB 
j a | 


M J s ° e e 
Where ee is the maximum variation among various 


<ag ms. 

Fig. (4-3) shows the noise disturbance behavior in 
the four Zi 5's and their average <Z> for raw data and 
for the indicated cyclic operations. Data was selected 
from sections of MT records with high noise contents. 
The result indicates that data between t ~ 10 to 100 gec 


are acceptable for the example shown. 
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Fig. (4-4) shows the result which is effected by 
the cyclic operation of Eqn(4-37). Data was taken from 
the MT record at the site-10 for both short and long 
period bands. Fig. (4-4.1) shows the average apparent 


resistivity cn (corresponding to <2 ey) of the raw data, 


y 
which have low predicted coherencies. Predicted coherency 
and smoothness of the curve are usually used as criteria 
for final data selection. If predicted coherency ¥ 0.9 
is taken to be the criteron, none of the data at this 
site will be accepted. 

To examine the effects of the cyclic processing, 
the low coherency of this data is ignored and a smooth 
curve is drawn through the data (Fig.4-4.1). Fig. (4-4.2) 
shows the average apparent resistivity after processing 
by the 6Evcliceoperation... Lf) Ch = 40/7) ~ 10% is taken to 
be the criteron, some of the data is now found to be 
acceptable, while other parts are still of no use. Using 
only this acceptable data gives an entirely different 
picture from the curve of Fig. (4-4.1). The results shown 


in figure (4-4.2) are also used in figure (5-5.10). 
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Fig. (4-3). Same as above for case when noise exists in 


both auto- and cross-power. 


Fig.(4-4.1). Average apparent resistivity om (corres- 


ponding to <Zgy?) OL) Law data tor both shorte(left) and 
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long (right) period bands. The numbers are the predicted 


coherency values for each indicated period interval. The 


curves are estimated by taking all the data into account. 


Fig. (4-4.2). Same data as above but after 10 cyclic 
operations. Acceptable data (CL 2108) are barred. The 
curves are established by using this acceptable data, 
and are also used for interpretation as shown in 
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CHAPTER V 
MT SOUNDING 


V-A. Geological Background 

The Black Hills are a pronounced surface feature 
which rise several thousand feet above the surrounding 
Great Plains of North America. The location is shown 
in Fig.(5-1). This island like mountain structure 
which is approximately 200 km long and 100 km wide 
extends North West from the South West corner of South 
Dakota into the North East corner of Wyoming. A con- 
siderable number of geological and geophysical studies 
have been carried out in this region and it is possible 
to give a general geological description. 

As shown in Fig. (5-1.2, 5-1.3, 5-1.4), the central 
Precambrian core of the Black Hills consists of a 
succession of highly folded schists, sedimentary in the 


main, but intruded by a succession of large and small 


granitic masses. Flanking this Algonkian mass are the 
upturned truncated edges of the sediments ranging in 
age from upper Cambrian to Tertiary. The dips of the 
sediments away from the centre are steeper on the East 
than on the West. Lidiak (1971) has mapped a meta- 
morphic belt under the sediments and trending somewhat 
East of the strike of,the Black Hills. This belt runs 


right through the central exposed metamorphic rocks of 


90 


{ me: ay AL 4 " 
H 
j i 
yer “7 } 
« 
t * 
™~ 
a . 
Fe v 
. 
4 , i 
1s b MG J 
H wv 7 
¥ ™ i a 
a r @ } 
‘ j ’ 7 
a 4 j eal : 
4 
x 
ae , 
i» ' 
_— 
i 
“ 4 i 
5 \ ¥ 
Le 
é 
F 
¢ 4 
a m 
= 2 Be iT cd eth 
. os 
“ 
a id t 
a * 
¥ 
7 
~ } \ 
- 
‘ ey 
7 * “ 
P «? .% " S 
Zi os 
bd +4 : 
' - : 
. " Sy 
3 4 : a . 
* @ 
‘ 
. a. " is 
co” 
‘ - « 
a i ’ 
=, =i 3 
es 


« i " 


” 


bs * % : Ke edt ie bs 
- ate DB 22 ‘ie 2 TUS 


u i on 
- “— 
Kina So 
= u 
I 
, 
— 
aM ® § 
: a 
‘ wd 
| 
f 
we x. 
aa Fs 
; oe 
uk i, 4 a 
lied 
} 
ji 
' il) 
ma 
o i 
d 
a 
A 
4 
5 = 


“ 
hide oF - 
( 
~ 
i deal 
i oJ Ed 
nim os 
PTA oa 
Me : 9 


<z 


erramde 
~ a LE ry 


ad 


e 


” 


OL 


the Black Hills, shown in Fig. (5-2). 

Gough and Camfield (1971) have found from their 
deep sounding studies a very pronounced induction 
anomaly which is in striking agreement with Lidiak's 
metamorphic belt, as shown in Fig. (5-3). They attribute 
their results to a suitable graphitic schist in the 
basement. Mathisrud and Sumner (1967) have found highly 
conductive graphitic schists in the lead mine district 
of the Hills, but these were highly localized. 

Three heat flow measurements made by Sass et al 
(1971), were either anomalously high or low. However 
these measurements were not of high quality and there- 
fore no reliable deductions can be drawn from them. 
Surface manifestations of thermal activity such as hot 
springs and elevated ground water temperatures exist in 
Many regions in and around the Black Hills. 

The age dating results quoted by Lidiak place 
Ghey blacks Hillis" orogeny at 1600) to. 1800 m-yeeagor 
Various silicic intrusives associated with volcanic 
activity are dated at 1450 m.y. ago. It should be 
noted that the extension of the metamorphic belt mapped 
by Lidiak was deduced on the basis of geophysical mea- 
surements such as gravity and magnetic and is reasonably 
consistent with drill hole results and the aeromagnetic 
survey reported by Zietz et al (1971) and also shown in 
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jer he a Co tow! lo Her Geological location of the Black Hills and 


sounding sites (@). Contour interval _2500-ft. Datum 
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Big.(o—-1.3). sthickness of Phanerozoic rocks: 


(Depth to Precambrian basement, in KM)-. 


Black Hills=uplitt (Precambrian crystalline rocks 


exposed). 
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Fig. (5-3). Contour map of Fourier amplitudes of 


east-west component of magnetic field at period 


48 min. Geomagnetic deep sounding result by Gough 
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Since electromagnetic measurements are affected 
by structure lying within several skin depths an aver- 
aging effect is anticipated and thus one would expect 
Lo measure a relatively, poorly. conducting, central 
megion.of the Black Hills characteristic of old, crys- 
talline material, surrounded by the relatively better 
conducting material characteristic of more recent | 
SSditients- wanemstrikingshigh conductivity, anomaly in 
the central region of the Black Hills will be discussed 


in the following sections. 


IV-B. Data Analysis and Computed Results 


Records were selected visually to be of good 
amplitude and free of spikes and steps. Six data sets were 
usually obtained at each station in both high and low 
mode. Details of this are shown in Table (5-1). In 
all cases a data set consists of 4096 digitized values 
with the high mode digitized at 40/sec and the low at 
1.25/sec. Following the techniques described in 
Chapter III, the Fourier coefficients were computed 
using a constant Q filter with Mb = 6 at the low 
frequency end and Mb = 300 at the high frequency end 
of “they spectrum. “Fig. (5-4) *shows?as pair of represen 
tative smoothed power density spectra for the E and H 


fields. The computed numerical results appear in the 


tables in Appendix B. 
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While it has been customary in the past to use 
only the criterion of high predicted coherency between 
the conjugate pairs of E and H as a basis for acceptance 
Ore rejeClion = Onedata, the cyclic cperation presented in 
Chapter IV permits the inclusion of a considerable 
amount of data with lower ntedicted coherencies in the 
initial stages of computation. The actual criterion 
adopted in the final stage of computation was that Pp. 
was derived from the average values of cn where the 
individual Ze lie within 24% of the average of the 
4 possible values. The results of this process can 
be seen in Figs. (4-2) and (4-3). By this procedure 
"good' data was obtained over the complete range of 
recorded frequencies, except for stations 3 and 9 where 
only a portion of the spectral range was acceptable. 

The apparent resistivities are shown in the 
Poeinet ot ecurectlons In P1lgse-(5— 10 GO) 5-5 .14) for 
each station. Real principal directions actually exist 
only for genuinely two-dimensional structures. These 
directions are given for each station as OF toe ka os. 
(Sat 51) Wey CHS 4 25) a 

In addition to the apparent resistivity, and its 
principal directions, the directional parameters OF 
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and oan are also used. These are the directions which 


produce a maximum coherency between He and a component 
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(L = long period band, S=short period band) 


Table 5-1. 
Site Period Band 
cl! Nahant L 
S 
Zz Custer L 
S 
5} Kahta L 
S 
4 Moon L 
S 
5 Smith L 
S 
6 Belle Fourche L 
S 
u Colony L 
S) 
8 Newell L 
S 
9 Wall L 
S 
10 Cheyenne River L 
S 
el Walker L 
S 
FIs. Ardmore L 
Ss) 
ks} Clareton L 
S 
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Fig. (5-4). Smoothed power spectral density. 
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Fig. (5-7.1). Contour map of maximum apparent resistivity 


(in 2-M) vs period along A-A' traverse. 
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of E and H respectively and are also shown in Figs. 
(5-6.1, 5-6.2) for each station. 

Pseudo-sections for the two traverses A-A' and 
B=-B" indicated in Fig. (5-1) are shown in Figs. (5-7.1) 
and (5-7.2), where apparent resistivity is plotted as 
a function of period. These are converted to MT 
resistivity cross sections in Figs. (5-8) and (5-9) 


and will be discussed in the next section. 


V-C. interpretation and Discussion 


(1) Discussions 


the results,/2.e. the. anisotropic character of 
the apparent resistivity and the dimensional indicators, 
show that in the area studied, the conductivity struc- 
ture is 3-dimensional. In the central region of the 
Hills, at sites 1-5, the two principal apparent resis- 
tivities diverge throughout most of the spectral range. 
The skew factor is large and also oe does not coincide 
with either hes or le as they would for a two dimen- 
ional earth. In the area external to the Hills, sites 
7-14, the apparent resistivity curves coincide for short 
periods and diverge for longer periods. This is inter- 
preted tominplyethatein) thisstegqionstheseart has 
homogeneous, isotropic, and layered near the surface 


and more complicated at depth. While the dimensional 


ie 


indicators are in disagreement with this interpretation 
it must be considered that these indicators are much 
more sensitive to the minor contributions from 3- 
dimensional structure. On the east, west, and south 
Sides around the Hills where sites 9-14 are located, 
the apparent resistivities indicate that below the 
conductive near surface layers, there is a relatively 
resistive layer, which in turn is underlain by a con- 
ductive third layer. In the area north of the Hills 
where sites 7 and 8 are located, it seems that the 
conductive third layer is not present. At the southern 
corner of the uplift at site 12, the principal axis of 
the apparent resistivity changes direction by 90° at 
tv10 sec as shown in Figs. (5-5.12) and (5-6.2); 

such a change only happens at this site. 


(2) Earth Model 


Taking site 2 as an example and applying the one- 
dimensional approach, an apparent resistivity curve 
aseshownean Hige(5-5.2) is used to construct a’ four 
layer earth with the parameters as given in Table (5-2). 
The first two layers would be acceptable due to the 
isotropic behavior of the apparent resistivities in 
the first decade of the period. The third layer has a 


very high conductivity and is located at approximately 


ie 7 : fs 
~ , ‘i 7 
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3-10 KM, corresponding to a period ranging of about 
1=500 sec. In this range the apparent resistivities 
appear anisotropic. Thus this anomaly may be placed 
at a distance of 3 KM from the measuring point but in 
an arbitrary direction in the half space. However, 
as shown in Fig. (5-7.2), for the same period range 


8 is closer to os rather than to of According 


fo) gs 
to the results of the two-dimensional model calcula- 
tions,6. coinciding with ces is consistent with the 
measuring site being on the conductive side of the 
strike of a resistivity anomaly. Thus, this low 
resistivity anomaly is placed below the measuring 
Site. This result is consistent with that of the 
group of stations in the anomalous region. 

The model parameters of layered earth for each 
site are calculated and listed in Table (5-2). The 


Maximum apparent resistivity Pu corresponds to the 


direction in which most of the power is concentrated 
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and thus Ou is considered to be more reliable for curve 


Bite ing mene msoULarlines: nwt OC oO) eco Te SCO CECE 
model results. 

Pseudo-depth, defined by z(w) = 6(w) are calcu- 
lated using the skin depth theory for each one- 
dimensional profile. The pseudo-depth represents the 
depth reached by the field with its amplitude reduced 


by a factor 1/e for a given frequency in a given 
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layered model. Contours of apparent resistivity and 
admittance vs pseudo-depth are referred to as MT 
resistivity cross sections. These are shown for the 
two given traverses in Figs. (5-8) and (5-9). 

While the MT resistivity cross sections are not 
adequate for detailed interpretation, they do give a 
qualitative picture which is useful in the total inter- 
pretation. 

Considering all the results, a final 3-dimensional 
resistivity section can be constructed. Figs. (5-10) 
show this construction for two traverses AA" and BB'. 
Figs. (5-11) show a plane view of the same resistivity 
distribution at horizontal planes at depths of 1-KM, 


and 5-KMjyand 20-KM, respectively. 


(3) Summary and Comments 


The cyclic operation developed in Chapter IV has 
proven very effective in extending the spectral range 
of data that can be utilized with a high confidence 
level. This has enabled a more detailed interpretation 
of the field data from the 1973 field season than had 
been heretofore considered possible. 

Tne the area surrounding the Blacky Hills upliett, 
low resistivity with some layering is found in rela- 
tively shallow depths underlain by material of high 


resistivity. Within the central region was found a 
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Table 5-2. 
Site R(Q-M) H (KM) 
Mahant 30-50 Aye) 
(1) 500 52:0 
5 Pept) 
2000 60% 
50 
Custer 30-50 aes) 
(2) 1000-3000 3. 
5 iE 
50 
Kahta 30-50 a5) 
(3) 1000-2000 8. 
50 4. 
50-400 
Moon 30-50 pe} 
(4) 2000-4000 re 
30 6. 
5:00 ZO. 
50 
Smith 30-50 age 
(5) 1000-3000 1S 


Belle 
Four. 


(6) 


Colony 
(7) 


Parameters of One-dimensional Model 


(R= Resistivity, H= Thickness) 


5000-10000 20. 
50.0 


10 48) 
2000-3000 .10. 
100 She 
200 

4 ate) 


2000-3000 20. 
100 


Site 


Newell 
(8) 


Wall 
(9) 


Chey. -R. 
(@mep) 


Walker 
(12) 


Ardmore 
(02) 


Clareton 
Gs) 


Morecraft 
(14) 


R (Q—M) 


4 
2000-3500 
200 

100 


Sie) 

5 
1000-2000 
100 
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100 


100 


50 KM 


allem, \@shasy) Contours of maximum apparent resistivity 
(upper, in 2-M) and maximum admittance (lower, in ar- 


bitrary units) vs z (pseudo-depth) along A-A' traverse. 
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piles) 


100 


100 


50 KM 


toe (590 Contours of maximum apparent resistivity 
(upper, in 2-M) and maximum admittance (lower, in ar- 


bitrary units) vs z (pseudo-depth) along BB' traverse. 
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Fig. (5-113) Resistivity (2-M) distribution in a 


horizontal plane at a depth of 20-KM from the surface. 
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thinner layer of low resistivity, a second layer with 
high resistivity and the closed low resistivity anomaly 
previously reported by Rankin and Reddy (1973), which is 
now interpreted to lie at a relatively shallow depth. 
While the low resistivity of the surface material is 
expected and usually attributed to the moisture and 
salt content of the recent sediments, the anomaly in 
the middle of highly resistive crystalline rocks is 
BiouemilErrculc Lomexpuain. since highly conducting 
material shields the deeper structure, the values 
assigned to the region below the anomaly are somewhat 
tentative but it appears that the low resistivity 
persists downward throughout the entire column of the 
crust. It is tempting to interpret this low resistive 
anomaly in terms of a thermal anomaly although other 
suggestions have been put forward. One other such 
suggestion by Camfield, Gough and Porath (1971) 
attributes their elongated induction anomaly to the 
existence of graphite schists in the basement. The 
aeQutteh las, that the anomaly is produced by a distri- 
bution of material with an abnormally high value of u 
cannot be overlooked and it is proposed to compute 
models embodying this concept. 

The extremely low values of apparent resistivity 
shown at stations 6, 7 and 8 are associated with peaks 


in the power density spectrum of the H field at 
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approximately 6.25 Hz. This is near the Schumann 
resonance peak reported at 8.6 Hz. This latter 
phenomenon is associated with waves propagating 
within ‘the earth-ionosphere cavity due-to: lightning 
strokes. If the peak described here is of the same 
nature as the Schumann resonance, the plane wave 
theory is not applicable and the results in this 
region of the spectrum may need modification. Work 


on this aspect of the problem is to be carried out. 
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APPENDIX A 


MODELLING CONSIDERATION 


Numerical modelling results have been published 
for numerous one and two dimension cases which have 
already been referred to. in this appendix several 


examples are presented schematically and discussed. 


I. One-Dimensional Model 


The typical feature of the one-dimensional solu- 
tions for a layered earth is the simple frequency 
dependence which uniquely reveals information on 
conductivity and thickness of the layers through the 
apparent resistivity and phase relations. [In Fig. 
(A-la), a pair of two-layer models and their apparent 
resistivity curves are shown. For sufficiently high 
frequencies, the corresponding skin depth is suffi- 
ciently small, so that the lower layer is not sensed 
and the apparent resistivity is asymptotic to Pi: For 
sufficiently low frequencies, the corresponding skin 
depth is large and the upper layer has little effect, 
so that the apparent resistivity approaches Po. The 
variation with frequency is smooth. A greater thick- 
ness of the upper layer, as in model-B compared to 


model-A, causes the apparent resistivity curve to be 
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displaced as shown. It is clear that the transition 
from Py CO P5 usa cunction of thesthickness Of ithe 
upper layer. 

Fig. (A-1lb) shows the result of a highly con- 
ductive and a highly resistive intervening layer. 
In both cases the apparent resistivity curves approach 
the same high and low frequency limits, while for 
intermediate frequencies the effect of the inter- 
vening layer is seen. The apparent resistivity never 
reaches the true values of P5 or ey as long as layer 2 
rs now very thick, but the effect of this» layer will 
always be present as long as it is not extremely thin. 
The extension to a multi-layer model can be easily 


visualized. 


II. Two-Dimensional Model 


Fig. (A-2) shows a solution for a vertical fault. 
Some notable features are as follows: 
a) The apparent resistivity is asymptotic to the 
appropriate resistivity value at large distances from 
the fault. It varies smoothly with position along a 
traverse crossing the fault in the case of Bi (E-field 
parallel to the strike), but varies discontinuously 
LOweche Hy case (E-field perpendicular to the strike), 


as shown in Fig. (A-2). 
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b) The vertical component of the magnetic field 

H, exists only for Bg with the maximum value 
appearing near the fault, decreasing to zero in both 
directions. The horizontal H-field perpendicular to 
the strike Hy also varies significantly near the fault 
for the ot case. Most of these effects can be 
understood by examining the current flow as shown in 
thes rELoure.. sFor ae and ED the current flows across 
the interface and since the current is continuous, no 
H, is induced. Since ein Sw ee, and EY, = OS .y1 Ay 
discontinuous across the fault with the condition 


E = E505. The apparent resistivity Pye which depends 


ea 
on Ee, changes abruptly across the fault and approaches 
the appropriate value at large distances on either side 
of the interface. 

Fig. (A-3) shows a pseudo-section for a fault. 
The horizontal scale is the distance along a traverse 
as used in Fig. (A-2), and the vertical scale is 
frequency (corresponding to depth) with the highest 
frequency at the top. The related data values are 
plotted beneath each site location for each frequency, 
and the values are contoured. 

Fig. (A-4) shows the effect of an overburden 
above the: fault seine result of ethesoverburden 1S eco 


smooth and attenuate the anomalous effects as shown; 


the degree of smoothing increases with the thickness 
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and conductivity of the overburden. 

Fig. (A-5) shows a conductive thin dike buried 
in the vertical position, and Fig. (A-6) shows a 
resistive dike. 'Thin' describes the condition where 
the dike is small compared to a wave length at a 
given frequency. The effect can also be understood 
Dyecxanining@chescurrene flow situation. For a con- 
ductive dike, the thin sheet-like conductive material 
immersed in a resistive basement does not affect the 
currents flowing perpendicular to the sheet; but it 
does play an essential role when the currents flow 
parallel to the strike. The resulting anomalous 
patterns for Pi) and He PSS WOW iiier ad Cm A— oo. a FOr 
the resistive dike, the situation is the converse; 
currents which flow across the dike are affected. The 


0 patterns are shown in Fig. (A-6). 


i 
Fig. (A-7) shows a two-dimensional model of a 


sloping contrast and the anomalous effects at a ty- 
pical frequency. Fig. (A-8) shows the anomalous 
effects for various frequencies at one location; the 


values are plotted as a function of frequency. 
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Fig. (A-1). a) Schematic response curves for 
two-layer models (upper). b) three-layer models 


(lower). 
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Fig. (A-2). Schematic pseudo-sections for a vertical 


Paulie. 


136 


O’ 
W 
4 
fy 
m Pu 
0 
4 
oF 
fe nt 
= 
m “s 
©) 
4 
o 0 0 
al Hy 
ey Hy 
S Max. 
Fig, (A-3). 2schematic pseudo-sections for 


a vertical -Eault. 


Loe 


Fig(a-4) Effect of overburden on a 
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APPENDIX B 


TABULATED ANALYSIS RESULTS 


PRD --- Period (sec). 

R12 --- Maximum apparent resistivity ({-M) for 0=0.. 
PH12 --- Phase angle of R12. 

R21 --- Minimum apparent resistivity (%-M) for 9=80. 
PH2L=-—— Phase angle of R21. 

REl --- Maximum admittance for 0=®0 7.. 

RH1 --- Maximum value of a for 9=0,,. 

ANG --- 0 5,from the geographical North. 

ANGE --- 96 

ANGH --- 954. 

SKW --- a, skew facter. 

BTA --- 85 (Eqn 2-54b). 

RAN --- 8] (Eqn 2-55a). 

AIS =-- 85 (Eqn 2-55b). 
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